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Abstract 

We investigate the observability of hydrogen molecules in absorption. The absorption efficiency of 
the hydrogen molecules become comparable with or larger than that of the dust grains in the metal- 
poor condition expected in the early Universe. If we can use bright infrared continuum sources behind 
the molecular gas clouds, the absorption measurement of the hydrogen molecules will be an important 
technique to explore the primordial gas clouds that are contracting into first-generation objects. 
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1. Introduction 

Molecular hydrogen is the predominant constituent of 
the dense gas in the Universe. One of the most important 
questions in current astrophysics is how the first genera- 
tion objects formed from the primordial gas that was not 
polluted by heavy elements (e.g. Omukai & Nishi 1998). 
Theoretical models should be confirmed by observations 
to establish a more realistic scenario of the first generation 
objects in the Universe. 

While in the local Universe, molecules containing heavy 
elements (e.g. CO, II2O) are good tracers of the amount of 
molecular hydrogen, they are expected to be significantly 
depleted in the early Universe, though some CO emission 
lines have been detected for high-z objects (e.g. Ohta et 
al. 1996). We need a technique for measuring molecular 
hydrogen directly. Petitjean et al. (2000) have reported 
a direct measurement of H2 molecules in the ultraviolet. 
Table |I describes some important transitions of hydrogen 
molecules. The transition probabilities of ionization and 
dissociation lines are so large that they are useful for de- 
tecting thin layers and small amounts of the molecular gas, 
but useless for detecting dense gas clouds. On the other 
hand, molecular hydrogen has well-known vibrational and 
rotational transitions in the infrared (IR) region. Their 
transition probabilities are very small because hydrogen 
molecule, a diatomic molecule of two identical nuclei, has 
no allowed dipole transitions but has allowed quadrupole 
transitions. 

These transitions of the hydrogen molecule can have 
both emission and absorption. Vib-rotational and rota- 
tional line emission from hydrogen molecules are useful 
tools to analyze dense (> 10 cm~'^) and hot (> 300 K) gas. 
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Expected emission line intensities of molecular hydrogen 
were calculated by Ciardi & Ferrara (2001). According 
to their result, direct measurement of the emission lines 
of molecular hydrogen is very difficult due to their weak- 
ness. On the other hand, if there is a strong IR continuum 
source behind or in the molecular gas cloud, absorption 
measurements of these transition lines may be possible. 
The absorption measurement has the great advantage that 
the column density of the absorber can be derived with 
fewer assumptions. 

However, we should note that the absorption by dust 
in the cloud is usually larger than the absorption by H2 
molecules. If the absorber is metal rich as in the ISM 
of the Galaxy, the extinction by dust causes the emer- 
gent continuum flux to decrease fast, thus making it dif- 
ficult to observe small absorption signals (Lacy et al. 
1994). Fortunately, metal abundance is expected to be 
significantly lower in earlier epochs of the Universe; hence 
the dust absorption is considerably less than that in the 
Galaxy. Therefore, hydrogen molecular lines can be de- 
tected in absorption against bright IR sources more easily 
in the early Universe than in the Galaxy. Such observation 
will be feasible with the advent of proposed space mis- 
sions for large IR telescope facilities, e.g. SPICA {HII/L2; 
Nakagawa et al. 2000), SPIRIT (Mather 2000). 

In this paper we investigate this possibility and con- 
sider what we can learn from such observations. The rest 
of the manuscript is as follows: In Section 2 we present 
the formulation of the intensity of absorption by hydrogen 
molecules. Our results are shown in Section 3. We discuss 
some relevant issues in Section 4. 

2. Calculation 

Assuming a uniform, cool gas cloud (fcTex <C hiy), the 
optical thickness of the line absorption, Tune is 
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Fig. 1. (a) Absorption line flux expected for a cloud with the Milky Way heavy clement abundance (Z = 1) in front of a 10 mJy 
source, (b) Same as (a) but the heavy element abundance is Z = 0.01. 



Table 1. Transitions of molecular hydrogen. 



Transition 


Wavelength 




A coefficients 


ionization and dissociation 


UV 




large 


vib-rotational transitions 


NIR (2 /im) 




small 


rotational transitions 


MIR (28, 17, •• 


■ Mm) 


small 


rotational transitions of HD 


FIR (112, 56, • 


• • Aim) 


intermediate 



rune ^ f [^-\ AuiA^-^ , (1) 
Stt \ g\ ) Av 

where the subscripts u and 1 indicate the upper and lower 
levels of a transition, t/u and gi are the degeneracy of each 
state respectively, A^i is the Einstein's A coefficient, N is 
the column density of the molecules in the lower state, and 
Av is the line width in units of velocity. Here we assume 
that almost all the molecules occupy the lowest energy 
state. As is well known, in the optically thin case, the 
column density can directly be derived from the equivalent 
width, which is nearly equal to the product of riinc and 
Av. 

The absorption line flux in the extinction free case, 
Atno.O' obtained by 

/,t.'^^,o-5Ai.(l-exp(-runc)) , (2) 

where Av is the line width in units of frequency and S 
is the continuum flux of the IR source behind the cloud. 
The line width and the source flux are assumed to be 
100 kms~^ and 10 mJy, respectively. We can scale the 
results of this calculation for other source fluxes, for ex- 
ample, based on the capability of the experimental set up. 
Assumed parameters are listed in Table ^. 

Next, we consider about the dust extinction which is 
denoted by 



W = 1.087 (^)(^)^ZiVH. (3) 

The extinction spectrum of Mathis (1990) is adopted for 
{A\/Av), and the extinction efficiency is assumed to be 
proportional to the relative heavy element abundance, and 
{Ay /iVn)© is the conversion factor from Ay to A'h for the 
local abundance listed in Table ||. 

Finally, we obtain the absorption line flux with extinc- 
tion, I^^l, as 

rabs rabs ^.^^^Z _ ^ 

^linc — ^lino,0'^^Pl~~dustJ 

= SAiy{l - exp(-Tiine))exp(-Tdust) • (4) 
3. Results 

The calculation was made for the five lines listed in 
Table 3. Parameters used here are from Turner et al. 
(1977). First, it is interesting to compare the optical 
thickness of the line absorption and the dust extinction. 
In case of Av = 100 kms^^ and the local heavy element 
abundance, the ratio of the line optical thickness to the 
dust optical thickness is 10~^ for vib-rotational lines in 
NIR, lO"*^ for pure rotational lines in MIR, and 10~^ for 
pure rotational lines of HD in the FIR. It means that the 
absorption measurement is certainly difficult against the 
large extinction in the absorbing cloud. However, this ra- 
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Table 2. Parameters assumed for the present calculation. 

Ortho : Para 
HD/H2 
Intrinsic Flux of Source 

Line width 
Detection Limit {5a) " 
Detectable Optical Thickness 
Dust Extinction Model 
a: Expected detection limit of the SPICA mission (Ueno et al. 2000). 



Table 3. Parameters of the Line Transitions. 





Transition 


Wavelength [/im] 


A Coefficient [s '-] 


A/ka [cm ^] for riinc ~ 0.01 


Ax/Av 


H2 


w = 1-0 S'(l) 


2.12 


3.47 X 10- '' 


4 X 10^^ 


0.11 


H2 


v = l-0 S{0) 


2.22 


2.53 X 10"^ 


8 X 10^3 


0.11 


H2 


S{1) 


17 


4.77 X 10-1° 


7 X 10^3 


0.020 


H2 


v = 0-0 S{0) 


28 


2.95 X 10"" 


3 X 10^4 


0.011 


HD 


v = 0-Q i?(0) 


112 


2.54 X 10-*^ 


2.5 X 10^4 


0.0011 



3 : 1 

10-5 
10 mjy 
100 kms-i 
5 X 10-21 Wm-2 
> 0.01 
Mathis (1990) 



tio increases with the inverse of the heavy element abun- 
dance. Therefore, in the lower heavy element abundance 
case, we can expect reasonably higher ratios. 

It is quite difficult to detect absorption lines whose op- 
tical thickness is less than 1 %. Therefore, as seen from 
Table 2, only those clouds whose column density is larger 
than lO^'i cm-2 can be detected in absorption. 

Figure la shows the absorption line flux expected for a 
cloud with the local heavy element abundance in front of a 
10 mJy source. The absorption fluxes are far smaller than 
the 5a expected sensitivity of the SPICA mission (Ueno et 
al. 2000). On the other hand. Figure lb shows the result 
for the case in which the heavy element abundance is 1 % 
of that of the local one. All five lines populate parameter 
space above or near the limit. 

Figure 2 shows the result of the same calculation, but 
for the II2 V = — 1 S{1) line for various values of the 
heavy element abundance. In case of Nn^ > lO^"' cm^^ 
and Z < 0.01, the absorption line can be detected. 

4. Discussion 

Possible Energy Source 

In this subsection we discuss possible energy sources 
at high redshifts. In our calculations, we have assumed 
an intrinsic continuum flux density of 10 mJy at the 
line wavelength. What sort of flux density can we ex- 
pect from high-z objects? To investigate this, we have 
plotted in Figure 3 the SED in the rest system of sev- 
eral high redshift objects for which observations are avail- 
able. The sources and references for them are as fol- 
lows: IRAS F10214-h4724 (z = 2.286; Rowan-Robinson 
et al. 1993; Barvainis et al. 1995), H1413-hll7 (z = 2.558; 
Barvainis et al. 1995), SMM 02399-0136 (z = 2.8; Ivison 
et al. 1998), SMM J 140 11 -I- 0252 (z = 2.55; Ivison et 
aL 2000), APM 08279-^5255 (z = 3.87; Lewis et al. 1988), 
and BR 1202 - 0725 (z = 4.69; Isaak et al. 1994). In this 



plot we have not corrected for the possible lens amplifica- 
tion for the first three objects and all observations have 
been transformed to a common redshift of 5 [qq = 0.1). 
The rest system SEDs and flux densities of F10214-f-4724 
and H1413-fll7 are almost identical and encompass three 
of the lines - 17, 28, and 112 /im. The interpolated flux 
densities at the three frequencies are 8, 16, and 9 mJy 
respectively. For the rest of the sources, we extrapolate 
from the observed 450 fiui flux density assuming the rest 
system SED to be the same as that of F10214 + 4724. 
From this exercise we find the range of flux densities to be 
1.6-19 mJy at 17 fim, 3-37 mJy at 28 /im, and 2-21 mJy 
at 112 fj,m. On the other hand, if we use the SED of 
Arp 220, the expected flux densities are about 38 times 
larger at 17 /im, 11 times at 28 /im, and about the same 
at 112 /im. For distant highly luminous sources, Arp 220 
SED may not be appropriate since Solomon et al. (1997) 
find that for a sample of ultraluminous IR galaxies out to 
a redshift of 0.3, the emission at 100 /im may be optically 
thick. Following their finding, if we use a 60 K blackbody 
spectrum as a template, the expected flux densities at 17 
and 28 /xm are a factor of 5 and 3 lower as compared with 
using Arp 220 SED. 

Thus, if we have objects at z = 5 at least as lumi- 
nous as our template sources, it is possible to observe 
absorption in H2 lines of 17, 28, and 112 /im respectively. 
The redshiftcd 112 /im line will be in the submm range 
and may not be observable by space missions of near fu- 
ture. However, the proposed large area and high sen- 
sitivity ground-based mm-submm interferomctric arrays 
like ALMA should be able to observe it. Unfortunately, 
it is very difficult to observe the absorption of 2.2 /tm 
lines since the expected flux density based on dust emis- 
sion spectrum is a factor of 10 or more lower as compared 
with that at 17 /im. However, if a bright non-thermal 
flat spectrum source is present, absorption measurement 
at 2.2 /im may be possible. 




Fig. 2. Result of the same calculation as in Fig. 1, but for the 
H2 V = — S{1) line in various values of the heavy element 
abundance. 

4-2. What can we learn from the absolution lines of pri- 
mordial clouds? 

First wc simply estimate the column density of a pro- 
togalactic hydrogen cloud. Here we assume that the size 
of the cloud i? is a few kpc and that the gas mass M 
is ~ 10^^ Mq. Such a large reservoir of molecular gas 
has been discovered at high redshift (Papadopoulos et al. 
2001). Such clouds have column density of 

A',..4X10-[C„-V(^)"(^).W 

where / is the mass fraction of the molecular clouds to 
the total gas mass, and M is the total gas mass of the 
protogalaxy. The final fate of such a large reservoir of gas 
may be the formation of a giant galaxy, whose typical core 
radius is a few kpc. If such a collapse really occurs in the 
formation epoch of galaxies, resulting column density is 
as large as the minimum column density (7 x 10^'^ cm^^) 
needed for detecting H2 line absorption. 

In a realistic situation, a primordial cloud may evolve 
dynamically on a free-fall timescale. The free-fall 
timescale is much shorter than the timescale of cosmo- 
logical structure evolution, e.g. the Hubble timescale. 
Therefore, observed properties are specific to the redshift 
at which the cloud absorption is measured. We will obtain 
redshift z, and velocity dispersion Av. The present struc- 
ture formation theory provides us the relation between 
the velocity dispersion and the corresponding dynamical 
mass of the cloud as a function of redshift, M^iriz) (cf. 
Padmanabhan 1993). Thus, we can constrain the struc- 
ture formation theory. The SPICA mission will provide us 
with the information on the massive objects 10^^ Mq) 
at z 5; 5. Wc still have to wait for more sensitive facil- 



Fig. 3. The spectral energy distributions of six representa- 
tive high-2 objects. Downward arrows represent the upper 
limits. 

ities for the observational approach to the "Population 
III" objects, whose physical properties are theoretically 
predicted by recent extensive investigations, because they 
are located at higher redshift and their typical mass is 
small (e.g., Nishi & Susa 1999). When we estimate the 
detectability of such objects by other facilities, the algo- 
rithm constructed in this paper can be applied straight- 
forwardly. 

The above estimate (equation 5) is based on the mean 
column density of a galaxy. If we have a chance to de- 
tect the final phase of a coUpsing gas in the line of sight 
of a quasar, we can constrain the physical state of such 
a coUpase. In other words, we will be able to measure 
the column density N^^, and excitation temperature Tox, 
which are interesting quantities in the context of the evolu- 
tion of an individual cloud. Recent theoretical work (e.g., 
Omukai & Nishi 1998; Omukai 2000) proposes the follow- 
ing scenario. The collapsing primordial gas must radiate 
its gravitational energy before the first generation objects 
arc born from it. However, the gas consisting of only 
hydrogen and helium atoms cannot cool efRciently below 
a few X 10'^ K because their constituents do not have radia- 
tive transitions corresponding a few x 10 K to a few x 10^ 
K. Therefore, a much more efficient cooling process must 
be working at this stage of the Universe. The most plau- 
sible cooling process is vib-rotational and pure rotational 
lines of molecular hydrogen. If the molecular hydrogen 
is effectively produced in the collapsing primordial gas at 
a few X 10^ K, the expected time scale of the gas contrac- 
tion is considerably reduced. Galli & Palla (1998) indi- 
cated that so called relic electrons would work there as 
effective catalyst to produce primordial molecular clouds. 
If we have a high-precision measurements and follow-up 
observations of the absorption cloud, we might be able to 
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pursue the dynamical evolution of the cloud and to com- 
pare the theoretical predictions. 
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